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STUDIES ON THE INTERFERENCE OF WINGS AND PROPELLER SLIPSTREAMS 

By 

Raaadas K. Prabhu 1 and Surendra N. Tlwari 2 
SUMMARY • 

The snail disturbance potential flow theory is applied to determine 
the lift of an airfoil in a nonuniform parallel stream. The given 
stream is replaced by an equivalent stream with a certain number of 
velocity discontinuities, and the influence of these discontinuities is 
obtained by the method of images. Next, this method is extended to the 
problem of an airfoil in a nonuniform stream of smooth velocity 
profile. This model allows perturbation velocity potential in a 
rotational undisturbed stream. A comparison of these results with 
numerical solutions of Euler equations indicates that, although 
approximate, the present method provides useful information about the 
interaction problem while avoiding the need to solve the Euler 
equations. 

The assumptions, of the classical lifting line theory applied to the 
wing-slipstream Interaction problem are scrutinized. One of the 
assumptions (uniform velocity in the slipstream) of the classical theory 
is dropped, and the governing equations are derived for the spanwise 
* lift distribution on a wing in a single axisymmetric slipstream. 
Spanwise lift and induced drag distributions are obtained for two 
typical cases, and the effects of nonuniformity in the slipstream 

Graduate Research Assistant, Deparment of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23508. 

z Eminent Professor, Department of Mechanical Engineering and Mechanics, 
Old Dominion University, Norfolk, Virginia 23508. 
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velocity profile are examined. 

The method of matched asymptotic expansions Is applied to the 
problem of a large aspect ratio swept wing In the slipstream of multiple 
overlapping propellers. The flow is assumed to be steady, Invlscld and 
Incompressible. It Is also assumed that the height of the slipstream Is 
of the order of the wing chord, and Its spanwlse extent is of the order 
of the wing span. Three different flow regions are Identified by 
employing different stretching transformations, and asymptotic 
expansions are Introduced using the chord-to-span ratio as the small 
expansion parameter. The details of the nonuniform flow In the 
slipstream enter Into the wing-sectional analysis. In the outer limit, 
the wing shrinks to a swept lifting line, and the slipstream reduces to 
a thin sheet of jet carrying the momentum gain from the propeller. The 
curvature of this jet sheet results in a pressure difference which is 
represented by a vortex sheet. The governing equations are solved by 
discretization. Comparison of the present results with the experimental 
data as well as other numerical solutions showed generally good 
agreement. 




Chapter 1 


INTRODUCTION 


The cost of aviation fuel has gone up substantially In the last 
decade and this Is expected to be an ever Increasing trend. Added to 
this Is the uncertainty regarding the supplies. These factors combined 
with the national concern over energy conservation have lead aircraft 
designers as well as the operators to give prime importance to fuel 
efficient propulsion systems for the future aircraft. Prior to the so- 
called oil crisis in 1973, the fuel cost was a relatively small fraction 
(about 25 percent) of the direct operating cost, and it was of less 
concern to the designers/operators. Today, however. It is claimed that 
this fraction has risen to about 60 percent [1, 2]* , and that it is the 
major part of the aircraft operating cost. 

. . . It was only propellers that, provided propulsive, force to aircraft 
before jet engines appeared on the scene. As flight speeds increased, 
the propellers posed serious problems of rapidly decaying propulsive 
efficiency and increasing noise and vibration levels. As a result, more 
powerful and efficient jet engines took over and dominated the scene, 
and propellers were neglected for many years. It is well known that the 
old technology propellers are the most efficient mode of propulsion up 
to a Mach number of about 0.6. The Interest generated in the propeller 


Numbers in brackets indicate references. 


technology since the oil crisis has lead to the development of the prop- 
fan, which Is claimed to operate at an aerodynamic efficiency of about 
80 percent at a Mach number of 0.8 [3]. 

With the prospect of the use of prop-fans on transport airplanes, 
there has been concern regarding the associated problems. The 
Interference of slipstreams with other parts of the airplane, in 
particular with the wing, is one of the major problem areas. With the 
flight speeds going up to H-0.8, compressibility effects can no longer 
be neglected. The flowfield behind the propeller is highly rotational 
and the effects of vertlcity In the slipstream cannot be ignored. 
Therefore, the problem in its complete form is quite complex. 
Consequently, considerable efforts (analytical, numerical and 
experimental) are being made to understand the rather complex flowfield 
associated with the wing-slipstream interference. 

The problem of determining the influence of the propeller 
slipstream on the wing lift has been studied quite extensively in the 
past, and a considerable amount of literature is available. Because of 
the highly nonuniform nature of the flow In the slipstream, the problem 
is essentially a nonlinear one; as a result, the work done during the 
1930's was based on approximate and semi-empirical methods [4, 5]. 
These methods provided satisfactory results In the speed-power range for 
which they were developed. Konlng [6] gave an analytical treatment for 
a wing in a propeller slipstream based on the lifting-line theory. It 
was assumed in this analysis that the Increment In velocity in the 
slipstream was small. Ferrari [7] developed what is generally referred 
to as the classical lifting-line theory for wings in slipstreams. The 
main assumptions in this theory were that (1) the slipstream was in the 



form of a circular cylindrical tube extending to Infinity both upstream 
and downstream, (2} the velocity In this tube was uniform (Uj), and (3) 
the relation between the lift and angle of attack for the wing sections 
was obtained by considering them to be in uniform flows with velocities 
Uj and U. for sections Inside and outside the slipstream tube, 
respectively. One of the drawbacks of this theory Is the third 
assumption. It is obvious that the lift produced by an airfoil would 
depend on the jet height. Ting and Liu [8] employed the method of 
Images and studied the lifting characteristics of thin airfoils in a 
nonunlform parallel streams. This method can be used to determine the 
lift of an airfoil In a uniform jet. Chow et al. [9] numerically 
Investigated the two-dimensional nonunlform flow past an airfoil by 
solving the Euler equations. These results demonstrated the fact that 
the lift of an airfoil depends not only on the jet height but also on 
the nonuniformity in the approach stream. Recognizing this fact, 
Klelnstein and Liu [10] made some improvements to the classical 
theory. The lift data for the wing sections within the slipstream was 
obtained by employing the methods of [8] and [9], and was used in the 
classical theory. However, the assumptions (1) and (2) of the classical 
theory were still retained while computing the downwash due to trailing 
vortices. These results demonstrated the effects of modifying the 

assumption (3) mentioned earlier. 

There were other attempts to improve upon the classical theory of 
Ferrari. Rethorst [11] employed the Weissinger approach and developed a 
lifting-surface theory. Wu and Talmadge [12], and Cumberbatch [13] 
extended the method of [11] to wings extending through multiple jets. 
Jameson [14] modeled wide slipstreams by rectangular and elliptic jets. 



and by using the standard Imaging technique developed a lifting-surface 
theory. In all these analyses one of the main assumptions was that the 
velocity In the jet (representing the slipstream) was uniform. 

Ting et al. [15] scrutinized the assumptions of the classical 

llftlng-line theory, and suggested a new approach to solve the problem. 
It was recognized that the height of the slipstream is of the order of 
the wing chord, and the spanwise spread of the combined multi -propeller 
slipstream is of the order of the wing span. Three different flow 

regions were identified by employing different stretching 
transformations. Asymptotic expansions were made by using the 
reciprocal of the wing aspect ratio (1/AR) as the small parameter. This 
analysis showed that the details of the nonuniform stream enter 

primarily into the local sectional anlysis; behind the wing, the 
slipstream acts like a thin jet sheet which supports a pressure 

difference across itself. By solving the governing integral equation, 
the wing lift distribution was determined. MaarsJngh [16] made an 
evaluation of this method by comparing these results with the data 
obtained from some specially designed experiments [17]. Some 
differences between the two sets of results were found, and these are 
suspected to be due to inaccuracies in the lift-curve slope data that 
was used in the computations of Maarsingh [16]. 

Ribner and Ellis [18] considered slipstreams of arbitrary cross 
section, and, instead of the standard Imaging technique, represented the 
slipstream boundary by vortex sheaths, and proceeded on the basis of 
lifting line theory, tan [19] developed a method based on the quasi- 
vortex-lattice method and a two-vortex-sheet representation of the 
slipstream. Both of these methods accounted for the rotation in the 


slipstream. 


Levinsky, et al. [20, 21] developed a large-tilt-angle lifting- 
surface theory applicable to tilt-wing and tilt-rotor V/STOL aircraft 
configurations. An actuator disk analysis for an inclined propeller was 
developed, and was combined with Weissinger lifting-surface theory for 
the wing at arbitrary wing angle of attack. Configurations with one, 
two, or four slipstreams were considered, and effects of slipstream 
swirl were Included in all but single slipstream case. Comparisons with 
experimental data showed that the theory predicts span loading 
reasonably well for small angles of attack and small propeller tilt 
angles. 

Kith the availability of high speed computers and efficient 
computational techniques, there has been a new trend In the approach to 
the solution of the wing-slipstream Interaction problem. Numerical 
techniques have been employed to solve the linearized potential flow 
equations, full potential equations, and Euler equations. 

Rizk [22] Investigated the propeller slipstream-wing Interference 
problem at transonic speeds. It was noted that a nearly uniform 
slipstream Interacting with a thin wing allows the perturbations to be 
potential although the undisturbed flow within the slipstream is 
rotational. The resulting potential flow (boundary value) problem for 
transonic flow was solved by a finite-difference scheme. In general, 
however, .the assumptions made in this work may not be fully valid; for 
example, the nonuniformity in the slipstream may be large enough to 
cause potential flow assumptions to yield erroneous results. 

Chandrasekaran and Bartlett [23] modified the Hess panel code to 


handle the effects of the propeller slipstream. The slipstream boundary 
was modeled by a system of ring vortices, and the effects of swirl In 
the slipstream was Included. A comparison of the results with 
experimental data showed some differences, which were attributed to 
viscous and compressibility effects as well as to uncertain! ties in the 
estimation of the flowfield behind the propeller. 

Naraln [24] and Samant et al. [25] made assumptions similar to 
those of Rlzk [22] regarding the slipstream, but did not assume that the 
perturbations were small. Instead, the problem was investigated on the 
basis of the full potential flow equations with a rigid boundary for the 
slipstream tube. These results compared reasonably well with available 
experimental data; however, it should be recognized that the advantage 
of using the accurate full potential equation is sacrificed by the 
approximation of irrotational flow. 

Whitfield and Jameson [26] solved the three-dimensional Euler 
equations coupled with the energy equation. The Euler equations had the 
force terms included to simulate the propeller effects. The .viscous 
effects were accounted for, although approximately, by coupling the 
three-dimensional Euler equations with the two-dimensional inverse 
integral boundary-layer equations. In spite of these sophistications, 
the spanwlse lift distribution obtained by this analysis failed to show 
good agreement consistently with the experimental data of Welge and 
Crowder [27]. This method, however, provides detailed information on 
velocity and other flew quantities in the entire computational domain. 
As may be expected, such an effort would require a large computer memory 
and a considerable amount of computing time. For example, the computer 
code developed in [26] required 900,000 words of memory and 341 seconds 
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of computing time on Cray-IS computer for a relatively coarse (96x16x16) 
grid. 

As a result of this literature survey, it is apparent that there is 
a need for 'further research in this field on the following topics: 

(1) Airfoil in nonuniform flow 

(2) Improvements to lifting line and lifting surface theories 

(3) Slipstream swirl and distortion effects 

(4) Compressibility and viscous effects 

(5) Swept wing - slipstream interference 

(6) Interference of slipstream with other parts of the airplane 

The main purpose of the present study is to investigate some of these 
topics. The following paragarphs describe briefly the work undertaken 
in this study. 

As noted earlier, the nonuniformity of the slipstream is not 
modeled properly in the classical lifting line theory for wing- 
slipstream interference. This has been Improved by using the sectional 
lift data obtained by either the linearized potential flow method for 
the wing section in an equivalent jet or the solution of two-dimensional 
Euler equations. Whereas the assumption of uniform flow for the 
computation of the sectional lift data In the classical theory is a 
drastic simplification, employing the Euler equations, although 
providing the necessary rigor, requires considerable computing effort. 
Approximating the actual nonuniform velocity profile in the slipstream 


by an equivalent uniform jet Is a good approximation. However, it Is 
found that better approximations are possible under the framework of the 
linearized potential flow theory. Chapter 2 includes a study of thin 
airfoils at small angles of attack in nonuniform parallel streams. 
These results are compared with numerical solutions of the Euler 
equations. 

It is assumed in the classical lifting-line theory for wing 
slipstream Interference that the slipstream is In the form of a uniform 
circular jet, for the computation of downwash due to the trailing 
vortices. The velocity distribution In the slipstream Is far from being 
uniform. However, this approximation (made to simplify the problem) was 
carried over In the subsequent developments. For example, Kleistein and 
Liu [10] retained this assumption in their modification of the classical 
theory, and so did Rethorst and his coworkers [11-13] in the development 
of the lifting-surface theories. In the present study, the assumption 
that the velocity in the slipstream is uniform is dropped, and the 
effects of the nonuniformities In the slipstream on the downwash of a 
large aspect ratio wing are studied. Chapter 3 Includes the lifting- 
line theory for large aspect ratio wings in slipstreams that can be 
represented as axisymmetric jets with smooth velocity profiles. 
Spanwise lift and induced drag distributions are computed for two 
typical cases. Results are compared with those obtained by other 
theories. 

The asymptotic method has been a powerful tool in the analysis of 
large aspect ratio wings [28]. The method Is simple, ard provides a 
better physical insight into the problem. The asymptotic method has 
also been used in an analysis of unswept wing-slipstream interference 
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CIS J. An asymptotic method for the analysis of the interference cn a 
large aspect ratio swept wing with multi-propeller slipstream is 
presented in Chap. 4. It is assumed as in £15] that the height of the 
slipstream is of the order of the wing chord, and that its width is of 
the order of the wing span. By employing different stretching 
transformations, three different regions are identified. Expansions are 
introduced in each region by using the ratio of the chord to the span as 
the small parameter. The details of the nonuniformities in the 

slipstream enter in the two-dimensional flow past wing sections normal 
to the lifting-line. The spanwise component of the velocity is shown 
not to affect the sectional lift data. For the outer solution, the wing 
planform reduces to a swept lifting-line, and the propeller slipstream 
behind the wing reduces to a thin sheet carrying the sectional momentum 
gained through the propellers. The curvature of this sheet results in a 
pressure difference across itself, which is represented by a vortex 
sheet. The governing equations are solved by a discretization 
procedure. Several examples are considered for which experimental 
results are available. Present results are compared with these 
experimental data as well as other numerical results. 


Chapter 2 


AIRFOIL IN NONUNIFORM PARALLEL STREAMS 

The study of aerodynamic characteristics of lifting surfaces in 
nonunlform flow is of considerable practical interest. Wing sections 
behind a propeller experiencing a jet-like velocity profile, and 
tallplane sections of a conventional airplane experiencing a wake-llxe 
velocity profile are two examples of such problems. These problems are 
complex and require simplifying assumptions. Even if the viscous and 
compressibility effects are neglected, the presence of vorticity in the 
approaching stream necessitates the solution of the Euler equations. 
Being nonlinear, the Euler equations require numerical treatment which 
has been done by several workers [9,29]. 

This nonlinear problem can be simplified considerably by replacing 
the- given nonunlform stream by an equivalent uniform jet. The advantage 
of this approach is the simplification of flow from rotational to 
Irrotational , with a finite number of surfaces of velocity 

discontinuity. Karman gave the basis for a linearized potential flow 
analysis of such problems [30], Glauert employed this method to solve 
the problem of an airfoil in the presence of a uniform jet [31]. The 
airfoil was replaced with a single vortex, and the airfoil lift was 
determined by computing the streamline curvature and increment in axial 
velocity. Ting and Liu [8] developed a method that is essentially an 
extension of Glauert's method, which could be used to compute the 
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chordwlse pressure distribution and the lift of a thin airfoil In a 
uniform jet. 

In this chapter the basis for linearized potential flow analysis 
for the problem of an airfoil near a surface of velocity discontinuity 
Is reviewed. The method of solution of the Integral equation of [8] is 
simplified. The analysis In then extended to cover the case of five 
streams with four surfaces of discontinuity. Next, the problem of an 
airfoil In a smooth velocity profile is treated by the linearized 
potential flow analysis. This problem is also solved by a more 

rigorous, although more expensive, method by solving the Euler equation 
using a modified version of the Euler code [29]. Results obtained for 
two examples by different methods are compared. 

2.1 Fundamental Basis for the Analysis 

The present linearized potential flow analysis Is based on the 
method due to Karman [30] for representing the flow past a body In the 
proximity of a surface of velocity discontinuity. Consider two parallel 
streams with velocities Ug and U^, the line AB being the undisturbed 
streamline separating the two streams (Fig. 2.1). Let a body be placed 
In the lower stream, and let (u 0 , Vg) and (uj_, v^ be the perturbation 
velocity components in the lower and the upper streams respectively. If 
the disturbed streamline separating the two streams makes an 

angle p with the undisturbed streamline, then 

tanp * v g/( U g + u g) “ V^/fUj+Vj) 

Assuming Ug, u^, Vg and Vj to be small compared to Ug and Uj, 
transferring this condition to the undisturbed surface, and retaining 



Fig. 2.1 A body near a surface of velocity discontinuity 
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only the first order terms, this relation can be simplified to 

VW U 1 «•». 

Also, the static pressure is assumed to be continuous across the 
surface, i.e., 

p Q - (1/2) pCUq - (U 0 +u 0 ) 2 - v 2 J - p x - (1/2) pCU 2 - (Uj+Uj) 2 - v 2 ] 

where p 0 and pj are the static pressures below and above the surface of 
discontinuity. Retaining only the first order terms, this relation 
reduces to 

u 0 U Q - u 1 U 1 (2.2) 

Equations (2.1) and (2.2) are the two necessary conditions that must be 
satisfied across the undisturbed surface of discontinuity, and form the 
basis for the analysis in this chapter. 

Glauert [31] considered an airfoil near a surface of discontinuity 
(Fig. 2.2). The airfoil was represented by a vortex of strength r. 
Since the problem is- linear, it was demonstrated [31] that application 
of Eqs. (2.1) and (2.2) lead to a flow in the upper stream as that due 
to a point vortex of strength (T+JOUj/Uq at the point P, and a flow in 
the lower stream as that due the vortex r at P together with its image 
of strength K at the point P‘, where 

k - r(u 2 - u 2 )/(u 2 + u 2 ) =• rp 

A logical extension of this approach is to replace the airfoil by a 
vortex distribution y{x), 0<x<c, Instead of a single vortex. Each of 
the vortex elements of this distribution y(x)dx forms images as 


(r+K)iyu 0 


Below AB Above AB 

Fig. 2.2 Image system for an airfoil near a surface of velocity 
discontinuity. 



described above. The downwash at the airfoil chora can be determined In 
terms of y(x) and its image strength hy(x) , and the unknown y(x) can 
be determined by satisfying the flow tangency condition on the airfoil 
mean camber line. This problem can be treated as a particular case of a 
more general problem of an airfoil in a jet of finite width, which is 
considered in the following section. 

2.2 Airfoil in a Jet of Finite Width 


A more Interesting problem Is the flow past an airfoil placed in a 
uniform jet of finite width. Consider the general case of th^ee 
parallel streams of velocities Uj, Ug and U_j with two surfaces of 
discontinuity AA and 0B separated by a distance h as shown in Fig. 2.3. 
Let an airfoil of chord c be pla >d in the middle stream at a distance a 
below the surface AA. In this case, the conditions (2.1) and (2.2) have 
to be satisfied at both the surfaces AA and BB. Ting and Liu 
represented the airfoil by a vortex distribution y(x)» 0<x<c instead of 
a single vortex [8]. The conditions (2.1) and (2.2) were applied 
repeatedly across the surfaces AA and 88, and an infinite set of image 
vortex distributions (Fig. 2.3) were obtained to describe the flow in 
the middle stream. The downwash v(x) at a point x on the y-axis due 
to y( x ) and all its Images is given by 


r (x) * ht f ^GFe) + 2 t J + 

6 0 IX u (x-sr + 4(jh+ar 


v(x-c) 


(x-O 2 + 4(jh+t>) 2 


+ ... . ]} yU) (2.3) 

(x-cr + 4(j+uV 


P “ (Ug - U 2 )/(Ug + U*) (2.4) 


where 


V* 




Part of the Vortex System 
for the middle stream 


Fig. 2.3 Image system for an airfoil in the middle of three 
; j parallel streams. 
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v ■ (Uq - U^J/OJq + U*) (2.5) 

The linearized flow tangency condition on the airfoil requires 

v(x) » UgC« - m(x)] on y«0 (2.6) 

where a. Is the angle of attack and m(x) Is the slope of the airfoil mean 
camber line. Equations (2.3) and (2.6) form an Integral equation for 
the unknown y(x). Note that In Eq. (2.3) the first term is the 
familiar singular term that appears in the classical thin airfoil 
theory; the other terms are not singular. 

The above equation was solved in [8] by employing a rather lengthy 
procedure. As the Integrand in Eq. (2.3) is no more singular than the 
one In the classical airfoil integral, all the methods of solving the 
classical Integral are applicable in the present case also. In 

particular, discretization of y(x) (s possible. Lan' s method of [32] 
discretizing y(x), and employing a cosine distribution for the vortex 
and control points is known to produce excellent results for the 
classical airfoil problem. Hence, this method is used In the present 
problem. 

As the first step towards the solution of Eq. (2.3), x and £ are 
replaced by 0 and $ using the following transformation: 

x * (1 “ cos 0) c/2 (2.7a) 

5 - (1 - cos$) c/2 (2.7b) 


Then, the Eqs. (2.3) and (2.6) together transform into the following: 


U 0 [a - m(0)J " 2 ! / { - ( C Qs ' ^co ~ e) ' + E C j — j 

0 ■ i* 0 lcos$-cose; K „ 0 (cos$~coser + 16(jh+a)Vc 2 
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v(cos$-cos0) 


2iiv(cos$-cos0) 


t]} y( ♦) 


(cos$-cos0) ♦ 16(jh+b) /c (coss-cosor + 16(j+ir h /c (2 g j 


Next, the vortex points and the control points are chosen as follows: 


- (2k-l)*/2N, k-l,2,...,N 


0^ ■ ix/N, 


i*l,2,...,N 


( 2 . 10 ) 


With this, Eq. (2.8) reduces to 


DoC.-oie,)] I 


+ 2 (uv) j [ 


pJCOS^-COSOj) 


vJcos^-cosOj) 


(cos^-cosO^r + lSfjh+aJVc^ (cos^-cose^r + 16( jh+br/c^ 


2pv(cos^-cos0.) 

+ jr 5 - 7 - 7 ]} y(4l) sinA . 1=1,2 N (2.11) 

(cos‘> k ”Cos0.j 1 + 16( j+lrhVc^ K K 

This is a set of linear simultaneous equations for the unknowns Y($ k )» 
k*l,2,...,N, and can be solved easily. The lift and pitching moment 
coefficients of the airfoil are given by 


Cj » ( «/N) 2^ y^) sin$ k 


( 2 . 12 ) 


c m « -(*/2N) 2 y(\) sin ( 1-cos ^) 


(2.13) 


Note that the series with index j in Eq. (2.11) converges fast, and 
hence can be truncated. It has been found sufficient to take N=* 15 and 
obtain good accuracy for the y[x) distribution. 


« 

< 





2.3 Airfoil in the Middle of Five Streams 

Consider the problem of an airfoil in the presence of an Infinite 
series of jets of the same width h. Let the velocity in the nth jet be 
denoted by U n , — <n<+*. Let an airfoil be located In the middle of the 
principle jet in which the velocity Is U Q . If the airfoil is 

represented by a vortex of circulation r, then the flow in any jet can 
be described by an infinite series of equispaced point vortices at the 
centers of each jet. The strength of these image vortices for the nth 
jet is denoted by K(n,s), -«<s<+«. In general, the following relation 
is true: 

K(n,n) - 0, n ¥ 0 (2.14) 

KCO.O) * T (2.15) 


By applying conditions (2.1) and (2.2) at the surface of discontinuity 
between the nth and (n+l)th jets, the following recurrence relation can 
be obtained for the strength of vortices: 


P n+ I K{n+l,n+s+l) * K(n,n+s+l)' K(n,n-s), 


“n “ <«K +1 >'<W+l> 


— <s<» 


(2.16) 

(2.17a) 


p n " /(1 'V (2 - 17b) 

Equation (2.16) can be solved in principle; but as pointed out by 
Glauert the solution is extremely complex [31], If, however, only five 
streams are considered as shown in Fig. 2.4, then the problem is 
simplified to some extent. When the conditions (2.1) and (2.2) are 
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Jet System 

Part of the Vortex System 
for the *1iridle Stream 

Fig. 2.4 Image system for an airfoil in the middle of five 
parallel streams. 
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applied to the four surfaces of discontinuity, four relations similar to 
Eq. (2.16) are obtained with n assuming the values 2,1,0 and -1 only. 
Oenoting K(2,s), K(l,s), K(0,s), K(-l,s) and K(-2,s) by I s , J s , r $ , K s 
and L s respectively, Eq. (2.16) can be written for the four different 
values of n as follows: 


S 2 ‘s.2 ' °s*2 ' °2 J l-s 

(2.18a) 

f l J s*l ' r sH ' \ r -s 

(2.18b) 

'o r S * Vb K - S -l 

(2.18c) 

B -I "s-1 ■ L s-1 • “-1 l - s -2 

(2.18b) 


where 


I s ■ 0 for s > 2, 


jj - o, r Q - r, k_j » o 


L ■ 0 for s < -2 
s 


The solution of the set of Eq. (2.18) can be obtained by substituting 

successively posi tive and negative values for s. However, with some 

algebraic manipulations, it is possible to obtain the following 

recurrence relations for r : 

s 

r s * “l r i-s + ff 2 r 3-s " “l^s^’ s > 2, (2.19a) 

r s “ a O r -l-s " “-1 r -3-s " ^) a -l r s+2’ S < ~ 2 ’ (2.19b) 

r 0 B T, Tj * c^r, and ^ =* -o^r 


Some of the values of computed using these relations are given 
below: 


h ' ‘ V-i r ' r - 2 

r 3 ■ ,0 2 4 - V^ ,r 

r -3 " - 0 -l4 )r 

r 4 " (< ^“l - ”o“2^ ' V-l l £ ,r ' r -4' etc - 

In general r n ■ a n r, -*<n<'» where a n 's are constants that depend only on 
Now, following the procedure adopted in the previous 
section the airfoil is replaced by a vortex distribution y(x), 0<x<c 
instead of a single vortex. Then the images will also be vortex 
distributions a n y{x). Note that it is not necessary to place the 
airfoil in the middle of the central Jet. The image system for the 
primary stream, when the airfoil is offset from the centerline is shown 
In Fig. 2.4. The downwash induced by the distribution y(x) and all its 
Images, at a point x on the y-axis is given by 

, C , ® , » (X-5) 

v{x) * vy + E 2 2 

2,5 0 tx « n— (x-?) 2 + 4(nh+a) 2 


• a,„(x-5) 

+ £ —5 j-j }y( l) 

n— (x-5) + 4n h 
n/0 


( 2 . 20 ) 


This Is the required expression for the downwash at a station x on the 
y-axis. The unknown y(x) in Eq. (2.20) is determined by employing the 
flow tangency condition at the mean camber line. 

As the first step in the solution, x and £ are transformed into 9 
and 4 respectively as In Eq. (2.7). With these, Eq. (2.20) is 
transformed into the following: ’ 





J-jL 


; , < ' .. « .ui'w;ji. ' 't»y,.’ K i VU», » w» )| J,- jym wjmwjtj 
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1 it j - a 2n+1 (cos<r-cose) 

tJ 0 ^ cos $“ cos6 > n— « (cos4rcose) 2 + 16(nh+a) 2 /c 2 

- . a, (cos 6-cos e) 

+ Z . .ii. — j } y( 6) sin 6 d6 

n— (cos6-cos9r + lEn^lr/c 
n^O 


(2.21) 


On discretizing y(6) and choosing 6^ and 0^ as the vortex points and 
control points respectively as in Eqs. (2.9) and (2.10), the Integral in 
Eq. (2.21) can be replaced by a finite sum 

N . , • a 2n+l^ cos ^k” cos V 


'<V ■•k £ « IcosvcosS,) * „l. 


(cos6 k -cos0j) 2 + 16(nh+a) 2 /c 2 


+ Z 


a, n (cos6. -cose.) 

* W-? ) t(U sin* . i”l ,2, . . . ,N (2.22) 

n— (cos6i,-cos0j 4 + 16n £ h 4 /c 
ntQ K 1 

The linearized flow tangency boundary condition on the y-axis at the 
control points requires 


v(e t ) • U 0 [«- mlty], 1-1.2 N (2.23) 

The, resulting set of simultaneous equations can be solved for the 
unknowns y(6j,), k»l,2,...,N. The lift and the pitching moment of the 
airfoil can be computed using the relations (2.12) and (2.13) 
respectively. 


Although the summation In Eq. (2.22) goes from -• to +«•, It Is 
sufficient to take only a few terms. This is because the Image vortices 
at greater distances from the airfoil are weaker and contribute very 
little to the downwash. It is found that about 15 terms In the 
summation are sufficient to provide better than 0.1 percent accurate 


results. 
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The two problems considered earlier, namely an airfoil near a 
surface of velocity discontinuity (Sec. 2.1), and an airfoil in a jet of 
finite width (Sec. 2.2), can be obtained as special cases of the present 
problem. For example, if U£ * Uj and Ug ■ U_j ■ U.g, then the problem 
considered in Sec. 2.1 is obtained, whereas if U 2 * Uj and U _2 * U_j, 
then the problem considered In Sec. 2.2 is obtained. Thus a single 
computer program (meant to solve the present problem) can be used to 
obtain the results for all the cases considered so far. 


2.4 Airfoil in a Stream of Smooth Velocity Profile 

So far the problem of an alfoll In a stream with a finite number 
(one, two, or four) of surfaces of velocity discontinuity was 
considered. This analysis can in principle be extended to the case of 
an undisturbed stream with a large number of velocity discontinuities. 
The solution of this problem would be, as pointed out earlier, rather 
complex. However, if the changes in the velocities in adjacent small- 
width streams are small, then an elegant solution can be obtained. This 
approach can. be used to solve for the pressure distribution on an 

airfoil in a stream of smooth nonuniform profile. 

Consider, as in the previous section, a large number of jets each 
of the same width h, the uniform velocity in the' nth jet being denoted 
by U n . Let the airfoil be placed on the axis of a jet in which the 
velocity is U Q (Fig. 2.5), and let the airfoil be represented by a 

vortex r . The discontinuity surfaces cause image vortices to be 

formed. The strength of these image vortices is governed by Eq. 
(2.16). It Is assumed that the variation of velocity in adjacent jets 

is small, i.e., (U n +i ■ U n ) * u n << U n , for all n. In this case the 
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expression for o n reduces to the following simple form: 


- u: 


“""777 


n+1 


2U u 
n n 


+ u: 


n+l 


2u: 


+ 2 U n U n 

n n 


I 


+ u n 


for u <<U 
n n 


(2.24) 


With j a n J«l, a first order solution can be obtained for Eq. (2.16). 
The resulting Image system for the primary jet is found to be following: 


r 0 “ r 


r 2n ‘ 0 


r 2n*l ’ V 


■ " a r 

n 


at y»0 

at y ■ (2n+l)h n > 0 
at y * (2n+l)h n < -1 


(2.25a) 

(2.25b) 

(2.25c) 

(2.25d) 


If the airfoil Is represented by a vortex distribution y(x) instead of a 
single vortex r , then the image system would be very similar (Fig. 
2.5). The downwash at a point x on the airfoil due to y(x) and all its 
images is given by 


1 c i ® a (x-5) 


- £ 


°n (x ‘ 5) 


n— 1 (x-E) 2 + (2n+l) 2 h 2 


} y( 5) d5 (2.26) 


Note that a n 's are defined in Eq. (2.24). For a given (thin) airfoil at 
a (small) angle of attack, the slope of the mean camber line is known, 
and as in the previous section, the unknown vortex distribution y(x) is 
determined by satisfying the linearized flow tangency boundary condition 




on the y-axis. Note that In deriving Eq. (2.26), It was assumed that 
the differences In velocities In adjacent jets are small. 

Analysis of an airfoil In an Infinite series of jets Is of little 
practical Interest. However, when the differences In velocities In 
adjacent small-width jets are small, the velocity profile may be 
considered as an approximation to a nonuniform smooth velocity 
profile. It Is possible to formally extend the present analysis to the 
case of an undisturbed stream of a smooth velocity profile by reducing 
the width (h) and correspondingly increasing the number of jets. For 
small h (*dy), u ■ (dU/dy)dy, the expression for a n reduces to 

: (2 - 2, > 


where U and dll/dy are measured at (2n-l)h/2 3 ndy - dy/2. The corres- 
ponding image is located at (2n-l)h * 2ndy - dy. In the limit as h 
tends to zero, the summations in the integrand in Eq. (2.26) are 
replaced by the corresponding integrals. With this, the downwash Eq. 
(2.26) can be rewritten as • ' 


,(x ' • Ti > 0 l T1RT 


_ e 1 dU (x-E) dy 
0* dy (x-E)* + 4y* 


♦J - fa = f - d 4-^(o « 

— U dy (x-EP + 4y^ 


(2.28) 


If the given velocity profile U(y) is even in y, and the airfoil is 
placed on the line of symmetry, then Eq. (2.26) reduces to the following 
simpler form: 

C 2(x-E) / c , 

0 U Qy (x-E) + 4y 


v (x ) 




Txhr ' z{ *~ z) L “0 dy 1 y(5) d5 (2,29) 








It can be shown that, for U and dU/dy of the order of unity, the 
linearized flow tangency boundary condition on y=0 is 


v(x) » U Q Ca - m(x)] 


(2.30) 


where Uq 3 u(y=0). 

For a given smooth velocity profile U(y) with U(y) »< 0, -«<y<», the 
integrals within the brackets in Eq. (2.28) or (2.29) can be evaluated 
using any standard technique and the unknown y(x) can be determined 
following the method described in Sec. 2.2. 

2.5 Results and Discussion 

The lifting characteristics of a flat plat airfoil in nonuniform 
flow are determined using the potential flow approach of Sec. 2.4. 
Also, the given nonuniform velocity profiles are replaced by equivalent 
stepped profiles, and the methods of Secs. 2.2 and 2.3 are applied for 
the computation of the airfoil lift. These results are compared with 
the numerical solution of the Euler equations. 

The lift of a cambered Joukowski airfoil in a nonuniform stream was 
studied by Chow, et al. [9] ana the results are available in [15]. In 
this example the velocity distribution in the undisturbed stream was 
assumed to have the following Gaussian profile: 


U(y)/U = 1 + a exp C-(y/d)^] 


(2.31) 


where “a" is the maximum excess velocity nondimensionalized 
using U^, and "d" is "a" measure of the spread of the nonuniformity. 
The value of "d" was chosen as c/1.81 and the value of "a" was varied 
from zero (uniform flow) to 1.0. The airfoil was placed on the 


centerline of this nonuniform stream. As the present method Is based on 
linear analysis, the lift-curve slope computed by this method does not 
depend on the airfoil camber. Hence for the present potential flow 
study, the airfoil is replaced by a flat plate. 

Figure 2.6 shows the lift-curve-slope of a flat plate 

(nondimensionalized using the corresponding value in uniform flow) in 
the nonuniform stream plotted against the parameter "a". The potential 
solution obtained for the flat plate in (1) equivalent uniform jet (Sec. 
2.2), (2) equivalent modified jet (Sec. 2.3), and (3) the given Gaussian 
velocity jet are also shown in the figure. Also shown in the figure are 
the results obtained from [15]. One of the observations from this 
figure is that for a given value of the parameter "a", the potential 
solution with the uniform jet approximation gives the smallest value of 
c. whereas the given nonuniform jet gives the largest value. The 

Euler solution gives a value in between these, and it appears that the 
modified jet approximation gives results closer to the Euler solution. 

The available Euler code [29] for uniform flow past an airfoil 
section was modified so that it could handle a nonuniform approach 
stream. Using this version of the code, the chordwise pressure 

distribution and the lift of a NACA 0012 airfoil in a jet of Gaussian 
velocity profile were computed. The value of the parameter "d" was 
chosen as 0.5c and the value of "a“ was varied. As in the previous 

example, potential solutions are obtained for a flat plate in (1) 
equivalent uniform jet, (2) equivalent modified jet, and (3) the given 
Gaussian velocity jet. All these results are shown in Fig. 2.7. The 
potential solutions display a behavior observed in the previous 
example. However, the Hft-curve-slopes -retained by the Euler equation 
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Ffg. 2.7 Lift curve slope of a flat plate In a nonuniform 
flow. 



32 


are comparatively smaller than In the previous example. 

The chordwlse pressure distribution (Ac p ) on the flat plate 
(obtained by the potential solution) and on the HACA 0012 airfoil 
(obtained by the Euler solution) at 5-degrees angle of attack are 
compared In Fig. 2.8. The nonuniform velocity profile is assumed to be 
Gaussian with a » 0.5 and d/c * 0.5. It may be observed from this 
figure that there is a reasonably good agreement between the two 
results. 

As the last example, lifting characteristics of a flat plate 
airfoil In a jet with a modified Gaussian profile are studied. The 
nonuniform velocity profile is assumed to be the following: 

U(y)/U. - 1 + aj exp {- (y 2 /d|) } 

- a 2 Cexp{-(y-d 3 ) 2 /d|> + exp {-(y+d 3 ) 2 /d 2 )] (2.32) 

with dj/c » 0.6, dg/c * 0.25, d 3 /c ° 0.155 and ag/a^ = 0.7. The factor 
a^ serves as the variable. 

. This .velocity profile,- sketched in Fig. 2.9, is typical of the 
axial velocity distribution generally observed behind actual 

propellers. The maximum velocity occurs not on the centerline but 
slightly away from it on either side. It is obvious that the 
approximation of a uniform jet would not be meaningful In the present 
case. Hence, this approximation is not considered. However, potential 
solutions were obtained for a flat plate in an equivalent modified jet 
and In the given nonuniform profile. The results are shown in Fig. 2.9 
along with the Euler solution for the NACA 0012 airfoil. It is observed 
from this figure that in this case the potential solution predicts 
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values which are slightly smaller than those obtained by the Euler 
solutions. 

An interesting observation can be made by comparing the results in 

Figs. 2.7 and 2.9. For the example shown in Fig. 2.7, the Euler 

solution gave values of c, smaller than the corresponding potential 

a 

solution, but this trend was reversed in the results for the example 
shown in Fig. 2.9. It was shown by the solution of Euler equations for 
an airfoil in nonuniform flow [9], that the stagnation streamline comes 
from a region slightly below the airfoil location. This implies that in 
the example considered in Fig. 2.7, the airfoil (placed on the 
centerline of a jet-like Gaussian velocity profile) would have a 

stagnation streamline with a total pressure slightly lower than the one 
corresponding to the centerline. Also, this streamline comes from a 
region having positive vorticity. In contrast, for the example 

considered in Fig. 2.9, the airfoil (placed on the centerline of a 

modified Gaussian profile) would experience a stagnation streamline 
having slightly higher total pressure than the centerline streamline, 
’and coming from a region of negative vorticity. These differences might 
have contributed to the reversal of the trends observed in Figs. 2.7 and 
2.9. 

Some differences are observed between the Euler solution and the 
potential solution. Among several reasons for this, the important one 
is that the potential solution does not account for the Interaction of 
the vorticity in the undisturbed stream with the airfoil. The effect of 
the airfoil thickness which is neglected in the potential solution could 
have contributed to the discrepancy. Yet another factor is the 

approximations in the numerical solutions of the Euler equations. For 



example, in the computer codes [9,29] the perturbations are assumed to 
be zero on the computational boundary. Considering these factors, the 
differences between the potential and the Euler solutions do not seem 
unreasonable. 

If the undisturbed stream has a smooth nonuniform velocity profile 
and an airfoil is placed in it, then the pressure distribution and the 
lift of the airfoil are to be determined by solving the Euler 
equations. Being nonlinear these equations are not amenable to 
analytical study and require numerical treatment. When a small 

perturbation approximation is introduced, the Euler equations or 
equivalently the vorticity transport equation reduces to the following 
linear partial differential equation (written in terms of perturbation 
velocity components): 

U(Uy - v x )x + U yy v =■ 0 (2.33) 

where U = U(y) is the undisturbed nonuniform velocity, and u(x,y) and 
v(x,y) are the perturbation velocity components assumed to be small 
compared to U(y). Subscripts x and y represent differentiation with 
respect to x and y, respectively. Since U(y) is assumed to be known, 
the above equation is a linear partial differential equation with 
variable coefficients. This equation together with the corresponding 
continuity equation has to be solved to determine the airfoil lift. The 
linearized potential flow solution obtained earlier obviously satisfies 
the equation (u y - v x ) =» 0. Hence, it is evident that this solution 
would also be a solution of the linerized Euler equation only if U yy is 
small. 

The concept of allowing the perturbation velocities to be potential 
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although the approach stream Is rotational Is not new. Rizk [20] 
employed this concept when studying the effects of a slipstream having 
nonuniform axial and rotational components of velocity on a wing. It 
was assumed that ' the undisturbed flow in the slipstream was nearly 
uniform and that the disturbances due to the wing were small. This lead 
to a result where the perturbation velocities could be described by a 
velocity potential. This approach, although approximate, allows 
obtaining basic information about the interaction problem while avoiding 
the need to solve the Euler equations. One of the assumptions of the 
classical lifting-line theory for wings in slipstreams (and some of the 
related works) is that the propeller slipstream is in the form of a 
uniform jet. This is a drastic simplification. Although the present 
analysis does not bring out the interaction of the vorticity in the 
stream with the airfoil, it is hoped that it would provide a much better 
approximation than the earlier ones to the actual problem. 
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A- MODIFIED LIFTING LINE THEORY FOR WING-SLIPSTREAM INTERFERENCE 

Consider the problem of determining the lift distribution on a 
large aspect ratio wing as Influenced by a single centrally located 

propeller slipstream. The classical theory [7] solves this problem 

making the following three assumptions In addition to those of the 

classical lifting line theory applied to large aspect ratio wings: 

(1) The propeller slipstream is confined within a stream tube of 
circular cross section in which the velocity Is uniform (Uj). 

(2) The relation between the sectional lift and angle of attack is 

the same as that of an airfoil In uniform flow (with 

velocities Uj and for wing sections inside and outside the 
slipstream respectively). 


-/ 
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(3) While computing the downwash, the stream tube representing the 
slipstream Is assumed to extend from upstream infinity to 
downstream Infinity. 

Figure 3.1 illustrates these assumptions in some detail. Although 
the assumption that the propeller sHostream is a stream tube of 
circular cross section is reasonable, the assumption of uniform velocity 
within the tube is not realistic. The slipstream behind a propeller has 
neither a uniform velocity distribution, nor a velocity discontinuity. 
The second assumption concerning the lift-curve-slope of the wing 
sections washed by the propeller stream is also not realistic. 
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These rather drastic simplifications of the. classical theory 
prompted several workers to study the problem In detail. Rethorst [11] 
studied this problem and developed a lifting-surface theory. However, 

It was assumed that the propeller slipstream was In the form of a 
uniform circular jet. Klelnstein and Liu [10] scrutinized the 
assumptions of the classical theory and Improved on one of the j 

assumptions. The lift-curve slope of wing sections in the slipstream 
was computed taking Into account the finite width of the slipstream. \ 

However, the assumption that the slipstream was In the form of a uniform ! 

jet was retained for the computations of the downwash. These results j 

brought out the effects of the assumption (3) above. ! 

In the present analysis the slipstream is assumed to be In the form | 

of an axisymmetrlc jet with a smooth velocity profile and without a 
distinct boundary. The relation between the sectional lift and the 
angle of attack is obtained from a local two-dimensional analysis. For I 

the purpose of computing the downwash due to the trailing vortices, the i 

slipstream is assumed to extend from far upstream to. far downstream. ! 

With these assumptions, the governing equations are derived for the lift 
distribution on a wing In slipstream. 
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3.1 Governing Equations for a Wing In a Single Slipstream 

The classical lifting-line theory given by Ferrari [7] for the wing 
slipstream Interaction problem Is an extension of Prandtl's lifting-line 
theory for large aspect ratio wings. The equation governing the 
spanwise distribution of circulation r(y) is 

T(y) ■ (1/2) U c(y) c^ (y) [a(y) - w(y)/U] (3.1) 





where c(y) is the wing chord, c. (y) Is the lift curve slope and a(y) Is 
the angle of attack; also U a Uj for |y|<R, i.e., for wing sections inside 
the slipstream tube and for (y|>R, i.e., for wing sections outside, 
R being the the slipstream tube radius. The downwash w(y) is given by 
the relation 


W (v) . 1 {/ S 

Wty Tii " s (y-n) 


1 r * dr(n) . . f R dr(n) _ rf " R . r s > dr(n) , 

« l, -s ^ ' 2 U ^ u , V 

I » _ 


|y|>R 


where Cj * (p 2 -l)/(p 2 +l), e 2 = (;*-l) 2 /( p 2 +l) and Uj/U^. It: ma y be 
recalled that in deriving these relations, the three assumptions 
mentioned earlier have been made. Further, whan the slipstream is 

absent, i.e., Uj * , the factors and e 2 become zero, and Eqs. 

(3.1) and (3.2) reduce to those of Prandtl's lifting-line theory. 


"If the jet representing the slipstream has a small excess velocity, 
i.e., Uj-U - « u«U a> , then terms of the order of (u/U J in and e 2 
may be neglected. In this case e^» u/U^ and e 2 » 0; as a result, Eqs. 
(3.2a) and (3.2b) are simplified to 


w(y) 



dr(n) 

(y-h) 



f R -JfiaL) 

-R (y-R^/n) 
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|y|<R (3.3a) 


1 r [ dr(n) 
T* * _ s (y-n) 


-f (/"“*/ )-S«2L). 

® -s R (y-Rvn) 


y|>R, 


(3.3b) 


Now consider a high aspect ratio wing with the propeller slipstream 
going past it symmetrically, as shown in Fig. 3.2. Let the undisturbed 






velocity distribution be given by U(y,z) » U(r) * U [1+F(r 2 )] where r 2 * 

OD 

y 2 + z 2 . Outside the slipstream (i.e., for r > R), F{r 2 ) » 0 and U(y,z) 

• V 

The flowfield inside the slipstream tube is rotational. 
Nevertheless, it is assumed that the perturbation velocity field due to 
the wing is irrotational . This concept of potential disturbances in a 
rotational background flow was employed by Rizk while considering the 
wing-slipstream interaction problem [223 . The effects of the swirl and 
compressibility were included and the resulting equations were solved by 
a numerical technique. 

Under the ' assumption that the disturbances are potential, the 
lifting-line theory is applicable and the wing is replaced by a lifting 
line. The circulation r(y) at a station y on the lifting line is given 
by 

r(y) = (1/2) U(y ,0) c(y) c A (y) [a(y) - w(y)/U(y,0)] (3.4) 

a 

The lift curve slope c^ (y) is determined by considering the 

a 

airfoil section in a stream of uniform velocity U for sections outside 
the slipstream, and by considering the airfoil section in a stream 
having a nonuniform velocity profile at the corresponding spanwise 
station for wing sections within the slipstream. In the present case, 
the wing section at the spanwise station y * yj would be in a stream of 
velocity U(y=y^, z). The lift curve slope for the wing section in this 
nonuniform stream is obtained by a two-dimensional analysis. This can 
be accomplished by solving the Euler equations, which require 

considerable computational efforts. A simpler, although approximate, 
method is the linearized potential flow method described in Chap. 2. 



For the sake of consistency as well as simplicity the linearized 
potential flow results are used here. 

Before proceeding to determine the downwash w(y) in the present 
case, it is useful to recall the results in the classical setting, where 
the velocity within the slipstream tube (a circular section of radius R) 
is constant. Consider a vortex (representing the wing trailing vortex) 
of strength y located at a distance ti from the center 0 of the circle 
representing the slipstream tube (Fig. 3.3). First, consider the case 
where | n j < R. By applying the interface conditions of continuity of 
pressure and streamline slope across the surface of the tube it can be 
shown [6], that the flew within the circle is described by a vortex of 
strength y at n together with its refracted image of strength (e^y) at 
the inverse point R*/n ; whereas the flow outside the circle is 
described by a vortex of strength (l-^)y at n along with an additional 
vortex of strength (p e^y) at the center of the circle. Similarly, for 
the case where the vortex is located outside the slipstream boundary 
(i.e., | rj| > R), it can be shown that the flow within the circle is 
described by a vortex of strength (1 - e^y located at tj. and the flow 

outside the circle is described by the vortex y along with its refracted 

p 

image of strength (e^y) located at the inverse point R /n, and another 
vortex of strength (c^y) at the center of the circle. These results are 
illustrated in Fig. 3.3. 

Now consider a propeller stream with a smooth axisymmetric velocity 
profile. For the purpose of analysis let this stream be divided into a 
large number of stepped, concentric annular cylinders of width £r (Fig. 
3.4). Let the axial velocities in the adjacent annular jets with the 
interface at a radial station r be U and U + u. Consider a vortex of 
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strength y located at Q (OQ ■ tj ). It Is easy to see that the 
difference u In the velocities of adjacent jets at the radius r results 
In an Image system as described in the previous paragraph. First 
consider the case where |n| < r. It can be shown that the flow In the 
region |y| < r Is described by the vortex at Q with its refracted image 
of strength (Cj y) at the Inverse point T (OT « r*/n ), whereas the flow 
in |y| > r is described by the only vortex at Q. Next, consider the 
case where Jn| < r. It can be shown that the flow in the region |y| < 
r Is described by the vortex at Q along with its refracted Image of 
strength (-e^y) at the Inverse point T. These results are illustrated 
in Fig. 3.4. 


The downwash at the spanwise station P (OP » y) due to the vortex 
of strength y located at Q (OQ ° n ) and Its image (whenever 
applicable) resulting from the surface of velocity discontinuity at the 
radius r is given by 

Aw(y,rj) = - {-i- + 2 ^"} jy|<r (3.5a) 

1 y-r / tj 

|y|>r (3.5b) 

« 

for the region fn|<r, and 


Aw(y,n) 


JL /JLl 

4 it 'y-n 

|y| <r 

(3.6a) 

^ 1 y-r 2 /r, 1 

|y|>r 

(3.6b) 


for the region |n|>r. In the limit as Ar tends to zero. 





u - u(p + Ar) - U(r) « (dU/dr) dr « U' dr 


so that 


e x - -u/U - -(U'/U) dr 


Letting r to vary from zero to R the following expression is obtained 
for the downwash at y due to a trailing vortex of strength y located 
at n : 


• - fc & ♦ [/'" l - j , 1 S'-TT 

0 |y| y-r /n 


} . |y|>h| (3.8a) 


0 |n| y-r /n 


} . |y |< j r»/ (3.8b) 


If r(y) is the unknown circulation distribution along the lifting line, 
then y * -(dr(n)/dn) dr). Using this in Eqs. (3.8a) and (3.8b), and 

integrating from tip to tip the following expression is obtained for the 
downwash w(y) at the spanwise station y due to the trailing vortices 
resulting from the distribution r(y) as influenced by the axisymmetric 


w(y) * / My,n) dn 


■ [/’ - {/' y ♦ n l(/f 1 - / . dr(n) 

-s -s y o (nj u y-r/n 

+ f. y {(/^ - L JS'-t—} <*r(n)] C 

y 0 |y| u y-r v q 

This equation along with the relation 


T(y) - (1/2) U(y ,0) c(y) c, (y) C«(y) - w(y)/(U(y,0)] 


(3.10) 
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form the required integro-differentlal equation for the unknown r(y). 
For a given wing, c(y) and a(y) are known; in addition, the velocity 
distribution in the propeller stream U * U(y,z) is assumed known. The 
sectional lift-curve slope can be determined by the method described in 
Chap. 2. With these informations, Eqs. (3.9) and (3.10) can be solved 
for the unknown r(y). It may be observed at this stage that the effect 
of the nonuniformity in the slipstream is twofold; it modifies the 
sectional lift-curve slope as well as the downwash distribution. 

3.2 Method of Solution 

A simple method of solving the Eqs. (3.9) and (3.10) for the 
unknown spanwise distribution r(y) is to assume it to be piecewise 
constant. This distribution results In a finite number (say N) of 
trailing vortices. The contribution to downwash from each of these 
trailing vortices can be computed easily using Eq. (3.8l. A summation 
of these individual contributions over all the N trailing vortices gives 
the required downwash w. 

First, it is convenient to transform y and n into angular 
coordinates 9 and $ by using the following transformations: 

y » s cose (3.11a) 

q ■ s cos$ (3.11b) 

Next, the trailing vortices are placed at the following N spanwise 
locations: 

$ k ■ (2k-l)*/2N k*l,2,...,N 


(3.12a) 



The strength of the trailing vortices at these stations are denoted 
by y(k) 4ns U m . The control points are chosen at the following N 
spanwise locations: 

0 1 - In/N 1-1,2,. ...N (3.12b) 

If l\| denote the circulation at the spanwise station y^ ■ s cosOj, then 

1 

' r(1) « 4« U E y(k), 1-1,2, ...,N (3.13) 

k-1 

This expression is used in the discretized version of Eq. (3.10) which 
is rewritten as 
1 

4ns U E Y(k) - (1/2) U(1) c(1) C. (1) [a(i) -w(i)/U(1)], 
k-1 a 

1-1,2, ;..,N (3.14) 

where U(i), c(i), c. (1) and w(i) are the velocity in the slipstream, 
a 

the wing chord, the lift curve slope and the downwash at the’ spanwise 
station y.. - s cos0 <t respectively. The computation of downwash w(i) is 
simplified considerably as there are only a finite number (N) of 
trailing vortices. ' The contribution to downwash from each of the 
trailing vortices is given by Eq. (3.8). Hence, the downwash w^d) at 
the control point i due to the N trailing vortices (together with their 
images) from one side of the wing centerline is obtained by summing the 
individual contributions. The result is as follows: 


w^i) 
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for { cos |< | cos 0^ | (3.15a) 





There Is a similar contribution to tfownwash from the trailing vortices 
on the other half of the wing, so that the total downwash w(i) at the 
control point i is given by 

N 

w(i) = u a Z '/{ k ) [G(i,k;U) - G(i,-k;U)], 
k=l 

i=l,2,...,N (3.171 
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Upon using this expression for w{i) In Eq. (3.14) the following set of 
simultaneous equations Is obtained for the unknown y(k) , k=*l,2,...,N: 

i 

Z Y(k) [l + n(i ) {G(i,k;U) - G(1»-k;U)}] 
k*l 

N 

+ z y(Ic)p(1) {G(i,k;U) - G(i ,-k;U) } - n(i) a(i) U(i)/li (3.18) 

k*i+l 

where p(i) ■ c(i) (i ) /8ns (3.19) 

a 

For a given velocity distribution U(r), the integrals in the function 

G(i,k;U) can be evaluated using any standard integration method. The 

lift curve slope c, (i) of the wing section is computed making a two- 
,l a 

dimensional analysis for the nonuniform flow past the wing section at 
the spanwise station s cosGj. The linearized potential flow method 
described in Chap. 2 is used for this purpose. With this information 
all the coefficients in the set of simultaneous equations (3.18) can be 
determined and the equations can be solved for the unknowns , y(k), 
k=l,2,...,N. Finally the circulation r(i), lift, and induced drag 
distribution, and other quantities are computed. 

3.3 Results and Discussion 

As the first example, a rectangular wing of aspect ratio 6.0 is 
chosen. The velocity distribution in the slipstream is assumed to have 
the following Gaussian profile: 

U(y,z) « Ujl+a exp{- (y 2 +z 2 )/d 2 )] (3.20) 

with a 3 0.5 and d/s ■ 0.3. The spanwise lift distribution on the wing 
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with this slipstream Is shown In Fig. 3.5, along with the lift 
distribution for the wing In the uniform flow for comparison. The 
figure Includes yet another lift distribution on the wing with an 
equivalent jet of uniform velocity. The velocity and the diameter of 
this equivalent uniform jet are obtained by equating the mass and 
momentum in the given slipstream to the corresponding values in the 
uniform jet. The wing In the given nonuniform slipstream analyzed by 
the present method produces a C L * 0.58 whereas for the equivalent 
uniform jet the C|_ a 0.57. Even though the two total lift coefficients 
are nearly the same, the lift distributions are very different. The 

present method does bring out the effect of nonuniformity of the 
velocity distribution in the slipstream on the lift distribution. The 
Induced drag distribution Illustrated in Fig. 3.6 also demonstrates this 
fact. 

As the second example, a trapezoidal wing of aspect ratio 6.67 and 
taper ratio 0.5 is chosen. The velocity distribution in the slipstream 
is. assumed as . . . • - • • • • * 

U(y,z) = [l+aj exp{- (y 2 + z 2 )/d 2 } 

-a 2 exp {- (y 2 + z 2 )/d| } ] (3,21) 

In this modified Gaussian profile, the maximum velocity occurs not on 
the axis but away from it. With 0<a2<l+a^, this distribution is a 
better approximation to the velocity distribution in the slipstream of a 
propeller. In the example chosen a^ * 0.6, ag = 0.75, dj/s = 0.3, and 
dg/s = 0.05. The spanwise lift distribution on the tapered wing with 
this slipstream is shown in Fig. 3.7. This figure clearly shows the 
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Chapter 4 

AN ASYMPTOTIC THEORY FOR THE INTERFERENCE 
OF SWEPT WINGS AND MULTIPLE SLIPSTREAMS 

The asymptotic method was employed by Van Dyke [28] for the 
analysis of large aspect ratio wings in uniform flow. This method was 
subsequently applied to wings with jet flaps by Kerney [33] and Takuda 
[34], Ting, et al. [15] applied this method to study the interference 
of unswept wings and multiple propeller slipstreams. As noted earlier, 
the asymptotic method Is simple and provides physical insight into the 
problem. The present chapter deals with an extension of Ting's method 
to swept wings, i.e., the application of the asymptotic method to the 
problem of interference of large aspect ratio swept wings and multiple 
propeller slipstreams. 

4.1 Mathematical Formulation 

Consider a large aspect ratio swept wing in a uniform flow (Fig. 
4.1). Propellers placed ahead of the wing produce slipstreams which 
flow past the wing. It is assumed that the wing geometry and the 
velocity distribution in the slipstream are known. The problem is to 
determine the lift distribution on the wing as influenced by the 
slipstream. 

In the present analysis, the flow is considered to be steady, 
incompressible and Inviscld. This amounts to a considerable 
simplification; however, the resulting problem is nonlinear due to the 
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CLANK NOT FILMED 




nonuniform flow In the slipstream, and the governing equations are the 
Euler equations. 

In the analysis that follows, two right-handed Cartesian coordinate 
systems are used. The (x,y,z) system has the x-axis In the direction of 
the free stream and the z-axls in the lift direction. The (x'.y'.z 1 ) 
system Is obtained by rotation of the first system through the quarter 
chord sweep angle A about the z-axis. It should be noted that the z'- 
axls is identical to the the z-axis and that the y'-axis lies along 
the quarter-chord line of the wing planform. 


In the absence of the wing, it Is assumed that the undisturbed 
velocity field can be written as follows: 


V(x/c Q ,y/s,z/c 0 ) « 1 U„ 

(ahead of the propellers) 

(4.1a) 

■ i U(y/s,z/c Q ) 

(behind the propellers) 

(4.1b) 


(for large Jz|) 

(4.1c) 


where i is the unit vector in the x-direction. This amounts to assuming 
that there Is a jump in the axial velocity across the propellers and 
that there is no swirl in the propeller stream. The velocity field 
around a propeller Is quite complex. There is inflow ahead of the 

propeller. There is a pressure jump and not a velocity jump across the 
propeller. The velocity field in the slipstream is periodic rather than 
steady because of the finite number of blades, and has the swirl 
component in addition to an increased axial velocity. In spite of this, 
the drastic idealization implied in Eq. (4.1) for the velocity field 
around the propellers is quite common, and springs from the fact that 
detailed treatment of the swirl and periodicity in the stream makes the 
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problem extremely complex. Attempts have been made by several workers 
to take Into account the effects of the swirl In a simplified manner by 
assuming that the swirl only introduces a change in the spanwise 
distribution of' the sectional angle of attack. In the absence of any 
superior approach, this simple method is adopted in the present study. 

The analysis starts by making the assumption that for large aspect 
ratio wings with multiple tractor propellers, the height (h) of the 
slipstream behind the propellers is of the order of the mid-chord c 0 , 
whereas the spanwise spread (21) of the combined propeller slipstreams 
Is of the order of the wing span b*2s. Thus there are two length scales 
Cq and b in the problem; hence, different stretching transformations are 
possible. Following the classical analysis for large aspect ratio 
wings, the outer region is introduced with the corresponding stretched 
variables x.y.z defined as follows: 

x ■ x/s, y = y/s and z ■ z/s (4.2) 

With x,y,z fixed and AR+°> (or e-*Q, where e=Cg/s), the wing shrinks to a 
line,, in .which. all the singularities that may be used to represent it 
are concentrated. The undisturbed velocity becomes 

+ + 

V «• i (ahead of the wing) (4.3a) 

+ - 

« i y a (behind the wing except on the cut z=0) (4.3b) 

The momentum gained by the stream in passing through the propeller is 
now contained in this thin sheet which acts very much like. a jet flap. 
The spanwise distribution of the momentum in this sheet is given by 
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J(y) “ pc 0 J U(U - Ujdz (4.4) 

All the details of the flow past the wing section are lost In this 
outer limiting process. To recover these details, the inner region Is 

A A A 

introduced with the inner variables x, y, and z defined by the 
following transformation: 

A A A 

X - x'/c 0 , y « y'/s and z ■ z'/c Q (4.5) 

AAA 

The following relations exist between the variables (x,y,z) and 
(x* ,y' ,z' ) : 

A A A A A 

x* * xcosA + ysinA, y' ■ ycosA - xsinA and z' * z (4.6) 

Now, the undisturbed velocity behind the propellers can be expressed as 

V(x/c 0 ,y/s,z/c 0 ) = i U(y/s,z/c 0 ) 

*► 

■ i UC (y ‘ cosA - x'sinA)/s,z'/c 0 ] 

, ■ * . A A A 

■ • . ‘ * • - * * f U[(ycosA - exsinA),z] 

A A A -U| A A 

* i [U(ycosA,z) - cxtanA -^(ycosA.z) +...] 

ay 

♦ A 

» i U(ycosA,z) for c << 1 (4.7) 

The undisturbed stream for the inner region, therefore, is given by (in 
the transformed coordinate system) 

V 3 [U(ycosA,z)cosA, U(ycosA,z)sinA, 0] (4.8) 
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The analysis of the flow In the propeller :tream behind the wing 
requires the introduction of another region (called the third region) 
with the corresponding variables x,y and z defined as follows: 

x » x/s, y ■ y/s and z * z/c Q (4.9) 

In this region, the undisturbed velocity is given by 
-► -► + 

V ■ i U(y/s,z/c Q ) = 1 U(y,z) (4.10) 

The three flow regions are shown schematically in Fig. 4.2. Equations 
(4. 3), (4. 8) and (4.10) represent the undisturbed velocity field for the 
three regions considered. The effect of the presence of the wing is to 
introduce disturbances in these flowfields. The disturbances in 
velocity and pressure depend on the small parameter c=Cq/s which, in 
turn, is related to the reciprocal of the wing aspect ratio. In the 
following sections, these disturbances are introduced in the flowfields 
of the three regions, and the governing equations are derived. The 
solution obtained in each region is matched with the others as described 
in Sec. 4.4, and a solution for the entire problem is obtained. 

4.2 The Inner Region 

The disturbances in velocity components and pressure for the inner 
region are expanded in power series of c as follows: 

u(x,y,z;e) = u^(x,y,z) + oi^(x,y,z) +... (4.11a) 

v(x,y,z;e) » v^(x,y,z) + ev^(x.y.z) +... (4.11b) 

A A A A A (0) * * * 1 j * * * 

w(x,y,z;e) =« w v (x,y,z) + ew' '(x,y,z) +... 


(4.11c) 







) 


ft. 

.1 

/ 


i < - ' 

],w- 


3 V<iy 


' / • ' 

r /. 


/ 




66 


A A A A A / Q \ AAA Af * \ AAA 

p(x,y,z;e) » p l ; (x,y,z) + ep v '(x,y,z) +... 


(4. lid) 


The undisturbed velocity components given by Eq. (4.8) and the 
disturbance quantities given by Eq. (4.11) are substituted in the 
continuity and the momentum equations. Upon equating the coefficients 
of like powers of e on either side of the equations, the following 
leading equations are obtained: 


;<«» * ;« 01 ■ o 

X z 


(4.12) 


•< / 


(UcosA + u (0) )ui 0) + w (0, (UcosA + u (0) ) A » -pi 0) /p (4.13a) 


z x 


(UcosA + u (0) )vi 0) + w^(UsinA + v^)„ ■ 0 
x z 


(UCOSA + U {0) )wl 0) + w* (0, w: (0) 
X z 


- pl 0) /p 


(4.13b) 


(4.13c) 


j. 


A (0) * * A 

If a stream function <i» (x,y,z) is defined such that 


u 




;i 0) '- UCOSA + u (0) 
z 


(4.14a) 


i ; 

'll 
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x 


(o) „ . ;(o) 


(4.14b) 


then Eq. (4.12) is automatically satisfied. Now, it should be noted 

A 

that Eqs. (4.13a) and (4.13c) do not involve the v-component of the 
velocity, and that by eliminating the pressure, they can be reduced to 
the two-dimensional vorticity transport equation 
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■*( 0 ) 

This Implies that u is a function of <J> only, or equivalently u Is 
A (0) 

constant along ^constant ( 1 .e . , along the projection of the 

* 

streamlines on the y=constant plane). Thus, the problem reduces to 
solving the Poisson equation 



,„((» + 7 ( 0 ). 


(4.18) 



This equation applies to the flowfield behind the propellers. In the 
plane Of the propel lers, the disturbances due to the wing section may be 
neglected. Consequently, the boundary condition at the propeller plane 
would be 

A 

4> (0> * / 2 .U dz (4.19) 

0 

The other boundary conditions are (1) far away from the airfoil u 
and w (0) tend to zero, (2) on the airfoil the normal velocity component 
Is zero, and (3) at the trailing edge of the airfoil the pressure is 
continuous (which is equivalent to the Kutta condition). Note that in 






this problem, y appears only as a parameter. 

The solution of Eq. (4.18), or the equivalent Eq. (4.12), (4.13a) 

and (4.13c), with the above boundary conditions can be obtained by a 

suitable numerical technique. The results would yield the values of the 
“( 0 ) “( 0 ) 

velocity components u and w everywhere in the field including 
the airfoil surface. 

A 

The y-momentum equation, Eq. (4.13b), may be rewritten as follows: 
[(UcosA + u (0) )-^x + J ( 0 , 4j(UsinA + v (0) ) = 0 


-^(UsInA + v (0) ) = 0 


(4.20) 


/ it 

' [I 
1 1 1 


where 4- = (UcosA + u^)-^r + w^A: . Note that - 5 - is the derivative 

as ax az 

along the projection of the streamline in the x-z plane. Equation 

A (0) 

(4.20) implies that (UsinA + v ) remains constant along the 

A A 

streamlines projected on the x-z planes. Since the airfoil is one 

“( 0 ) 

such streamline it is concluded that (UsinA + v ) is a constant on the 

.airfoil,. Let U* be the velocity far ahead of the. airfoil on the zero 

* A (0) 

streamline. Then, it is evident that U slnA = (UsinA + v ). Note 

A (0) 

that, even though (UsinA + v ) is a constant on the airfoil, both 
A (0) A 

'slnA and v being functions of z, vary along the airfoil surface. 
The pressure coefficient c p on the airfoil defined by 


c p - ' (p-pJ/(l/2)pUt 


(4.21) 


can be reduced to 


c p (x,z) = [U* 2 cos 2 A - (UcosA + u<°>) 2 - (w l0) ) 2 ]/U 2 


(4.22) 
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Hence the difference Ac. Is given by 

P 

Ac p (x) * ‘ c p f x » z u ) 

« - {(UcosA + u (0) ) 2 + (w (0) ) 2 }/U^ (4.23) 

where { } represents the difference between the lower and the upper 
surfaces of the airfoil. The sectional lift coefficient c^ Is obtained 

A 

by Integration of Ac p (x) along the chord. Thus, It Is seen that the 
lift obtained by the solution of Eq. (4.13) Is not affected by the 
spanwlse component of the velocity. In summary, the governing equations 
for the inner solution are the continuity and momentum Eqs. (4.12) and 
(4.13), with the appropriate boundary conditions. The solution of these 
would yield the lift of the airfoil which is independent of the spanwise 
velocity component. This lift is a function of the sectional angle of 
attack, and may be written as follows: 

A{ a,y) = (1/2) pU 2 c(y)c A (a,y) (4.24) 

* * A 

Notice that in this relation y appears only as a parameter. 

4.3 The Third Region 

In the analysis of the third region, i.e,, in the propeller stream 
far behind the wing, the stretched coordinates x,y and z defined by Eq. 
(4.9), and the undisturbed velocity defined by Eq. (4.10) are used. The 
perturbation velocity components and the pressure are expressed in power 
series of e as follows: 

u(x,y,z;e) = u^(x.y.z) + eu^(x,y,z) 




w(x,y,z;e) = ew {1) (x,y,z) + e 2 w (2) (x,y,£) + ... (4.25c) 


p(x,y, z;c) « cp (1) (x,y,z) + e 2 p {2) lx.y.z) + ... (4.25d) 

These perturbation components and the undisturbed velocity given by Eq. 
(4.10) are substituted in the continuity and momentum equations. The 
subsequent analysis would be identical to the analysis performed in [15] 
while studying the third region behind unswept wings. Therefore these 
details are not given here; but only the major conclusions from this 
analysis are listed below: 


(1) To the first order, there is no pressure discontinuity or change 
in streamlinq inclination across the jet sheet. 

(2) The momentum integral J(y) remains constant with respect to x. 


(3) To the second order, there is a pressure difference across the 
jet sheet given by 


Ap (2) » p (2) (x,y,0+) - p (2) (x,y,0-) =■ - eeL 1) (x,y,0)J(y) 


(4.26) 


3 w (1) /u 


Equation (4.26) is the familiar result of the thin jet approximation 
that states the pressure difference across the sheet is equal to the 
product of the momentum in the jet sheet and its curvature. These 
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The undisturbed velocity in the outer region is uniform (Eqs. 
(4.3a) and (4.3b)); hence, the fiowfield away from the wing and the 


; , \ propeller streams is irrotational. Therefore the disturbance velocity 

l '( ■ fie1d can be described by a velocity potential &(x,y,z;e) which 

* \- satisfies the Laplace equation 




+ 4*r) ■ 0 

dz* 


(4.27) 


Bernoulli's equation provides an expression for the pressure 
disturbances as 

p = U/2)pu£ci - |i + WUJ 2 ] < 4 .28) 

Equations (4.1) and (4.2) are the governing equations in the outer 

region.. As the first step towards the solution, $ is expanded in a 
power series of e as 

®(x,y,z,*e) = e®^(x,y,z) + e 2 <S^(x,y,z) + ... (4.29) 

where both ^ and in turn satisfy the Laplace Eq. (4.27). Next 
the velocity and pressure disturbances are also expanded in power series 
of c as follows: 



u(x,y,z;e) = eu ll) (x,y,z) + e 2 u~ (2) (x,y,z“) + ... 
v(x,y,z;e) = cv (1) (x,y,z) + e 2 v (2) (x,y,z) + ... 


(4.30a) 

(4.2GL) 




• •*•*•* *•/ 1 \ •••*•* «» « « 
w(x,y,z;e) ■ ew l ; (x,y,z) + eV* ; (x,y,z) + 

p(x,y,z;e) » ep (x.y.z) + c p l '(x,y,i' + 

Thus, 

u< 1} » 4>i 1 *(x,y,z)/s 
x 

v (1) - “ «i^(x,y,z)/s (4.31b) 

y 

w (1) ■ sl^tx.y.zl/s (4.31c) 

z 

p (1) = pU. «l l) (x,y,z)/s (4.31d) 

x 

“( 1 ) 

The unknown ' is a solution of the Laplace equation. It is obtained 
by distributing the singularities along the line x-m|y|=0, and on the 
plane z=0, x-m|y|>0, and by matching with the solutions of the inner and 
the third regions respectively. Note that m is the tangent of the wing 
quarter chord sweep-back angle. 

By an analysis of the inner region it was shown that the wing 

section at a spanwise station y produces a lift force A(a,y) (Eq. 

4.24). This implies that there is a vortex distribution along the 

lit line x-m|y|=0. The strength r(y) is related to the 

lift Jt(a,y) as follows: 

T(y) - *(a,y)/p 

= {l/2)U o c(y)c x (a,y) (4.32) 
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(4.30c) 

(4.30d) 

(4.31a) 
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The spanwise variation of r requires the presence of a trailing vortex 
sheet extending downstream from the lifting line. This is the well 
known vortex system of lifting-line theory. 

By a matching with the solution of the third region, it was shown 
that, to a first approximation, there Is no discontinuity in either the 
pressure or the streamline inclination across the jet sheet behind the 
wing. Therefore, the first order outer solution is not affected by 
the propeller stream behind the lifting line. Based on this first order 
analysis, the flow in the outer region is described by the potential due 
to r(y) and the associated trailing vortex system, and is given by 

® {1) = J- / { j zr( . J } — =* {1 + dn (4.33) 

-1 (y-n) 2 + 2 2 R 

where R = C(x-m| n| ) + (y-n) + z ]. Upon approaching the lifting 
line (i.e., for small x-m|y| and z), the inner limit of ' is 
obtained. This yields the downwash 6(x-m|y|-*0,y,0). By matching with 
the inner solution, the effective angle of attack at each spanwise. 
station is obtained as 

- -Ml — ~ 

o(y) * a g(y) “ e ® U 3 <*g(y) ■ 9 (4.34) 

The circulation r(y), given by Eq. (4.32), now becomes 

r(y) « (l/ajU.ctylc^a - 0; y) (4.35) 

Equations (4.34) and (4.35) are the required governing equations for the 
first order analysis. Note that in this analysis, the effect of the 
nonuniform flow in the slipstream enters only in the sectional analysis. 



The first order analysis discussed so far is valid 

for e*0, or AR+® . It was shown in the analysis of the third region 

that terms proportional to J(y) appear in the se.ond order analysis; 

•* 

hence it is necessary that J{y) be of the order of 1.0 for the first- 
order theory to be valid. It is well known that the first order theory 
for large aspect ratio wings in uniform flow gives satifactory results 
for wings of aspect ratio as low as 1.0 [33]. Hence, only the 

contribution of the momentum in the propeller stream will be considered 
in an extension of the present analysis. 

“ill - 

At the outset It is noted that <&' is an odd function In z. It 
was shown by an analysis of the third region that p^ (x,y,z) Is 

continuous across z=0. Hence It is concluded that 

p U) (x.y»z=0) =0 (4.36) 

In view of Eq. (4.31a) and (4.31d), it may also be stated that behind 
the wing 

u (1) (x,y,z=0) = 0 (4.37) 

A* M tm 

However, because of the vortlclty dr(y)/dy on the plane z=0 behind the 
wing, v ll) is discontinuous across the jet sheet and is given by 

v^(x,y,z=0+) 3 dr(y)/dy « - v^(x,y,z=0-) (4.38) 

1 ) 2 

Nevertheless, [v ] is continuous across the jet sheet. Furthermore, 
since 8^ Is continuous, is also continuous across the jet sheet. 

2 

Next, upon considering the 0(c) terms in the Bernoulli equation 
the following equation is obtained: 
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p (2) /p - - U„u {2) - t(v (1) ) 2 + (w {1) ) 2 ]/2 (4.39) 

Since [v ] and w are continuous, the pressure difference across 
the sheet is given simply by 

« <2 >. p' (21 (x".;,;.o*) - p !2l (x'.;.;.o-) 

■ -pU a> [u^ 2 ^(x,y,z=0+) - u^(x,y,z=0-)] (4.40) 

-(2) 

This discontinuity in u is equivalent to a vorticlty 
«• «• 

distribution y(x,y) on the jet sheet and, in view of Eq. (4.26), 

M •> 

y(x,y) can be expressed as 

y(x,y) = - e£i 1 *(x,y,0)J(y)/pU ci s « - (Ux.y.OJJfyJ/pU^s (4.41) 
x x 

At the lifting line, 8 is not zero (i.e., the jet sheet has a non-zero 
Inclination at the lifting line). This would require a vertical force 

- ~(i) 

of magnitude J{y)c8 at the lifting line. To account for this force, 
a circulation distribution Tj(y) is introduced along the lifting line. 
The magnitude of this circulation is given by 

Tjfy) = - J(y)e0 (1) /pU a , = - J(y) e/pU^ (4.42) 

Thus, in this limited second order analysis, referred to as the 
systematic analysis in [15], the vortex system consists of the 

■a M 

circulation T(y) along the lifting line, the circulation Tj(y) which 
accounts for the change in the inclination of the jet sheet on passing 
over the wing, and a distribution y(x,y) behind the wing, which accounts 
for the curvature of the jet sheet. This system of vortices produces a 
downwash distribution 9(x,y) which is obtained from the following 
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velocity potential: 
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«(.x,y,z) - 4- / ■ ~- ^ r l Tl) . u + 


411 -1 z” 2 ♦ (y-n) 2 


dn 


1 ^ 2Pl(T)) Y-mlnl 


-l z + (y-n) 


+ 1 / /" *li|a2> u + (x ^ } 

’* "l. t, l z + (y “ T,) / C(x-0 2 + (y-ri) 2 + z 2 ] 


} d5 dn (4.43) 


The downwash Is related to the unknown quantities by the following 
relations: 

T(y) » (l/ajU^cty )c . Cat (y) - 6(x,y)], x + m|y| (4.44) 


«* «■ «• 
Tj(y) = - vHyjefx.yJ/pU^, x m|y 


(4.45) 


y(x,y) = - J(y) 0-(x,y )/ pU ot s 
x 


(4.46) 


These are the governing equations for the unknowns r(y), Ij(y) and 
r(x,y). . 


4.5 Method of Solution 

The governing equations for the unknowns r(y), Tj(y) and y(x,y) 
are the Eqs. (4.44) - (4.46) together with the Eq. (4.43) for the 
downwash angle 9 required in the solution. The downwash angle computed 
by Eq. (4.43) on the lifting line at the wing centerline is infinite 
because of the discontinuity in the slope of the lifting line at the 
wing root. This is a familiar problem in the swept wing analysis. 



Rigorous studies of swept wings in uniform flow by Cheng, et si. [36,37] 
have shown that the induced velocity apprdaches infinity like the 
logarithm of the distance from the centerline. This behavior does not 
allow computation of the spanload at the wing centerline of 
symmetrically swept wings. Thurber [38] studied swept wings with curved 
centerlines having zero sweep at the wing centerline and zero tip chord 
(crescent wings) in uniform flow. This planform has limited application 
in the present context. 

The method developed by Lan [39] is employed here to compute the 
downwash from Eq. (4.39). This method starts with the Weissinger 

velocity potential, and, placing the lifting line on the quarter-chord 
line, computes the downwash at the three-quarter-chord line. A brief 
description of this method is presented here. 

Consider the first integral on the right-hand side of Eq. (4.43). 
This is the potential due to r(y) along the lifting line, i.e., ' 


„ _ 1 zr(n) „ „ (x-S) 

l . 4* _l z 2 * * S + (y-n) 2 R 


)dn 


(4.47) 


2 2 2 2 

where R = (x-5) + (y-rj) + z and tildes are dropped for 

convenience. The origin for the co-ordinate axes is assumed to be at 
the mid-root chord. Also, 


S ■ - Cq/ 4 + m| n 


(4.48) 


x = - c q /4 + c(y)/2 + m[y | 


(4.49) 


These are Introduced in Eq. (4.47), and the resulting expression is 
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expanded for small c(y). By retaining only first order terms in the 
expansion, and differentiating them with respect to z, the required 
expression for (evaluated at z=0) is obtained as 


Mx.y.oi - ij / -SiaL a * aiM * M 
1 4lt -1 (y-n) l Rj 1 


(4.50) 


where m' ■ iii-{Cq/ 4){1-X) is the tangent of the mid-chord sweep, and 
R 2 ■ (x-m'|nj) 2 + (y-n) 2 . The right hand side of Eq. (4.50) is 
Integrated by parts to arrive at the following expression for Gp 


(i i. .. m , If 0(h) dr(ri) j . 1 f ri tw vi. \ ii < A _ 

0i(x,y.°) - ^ dn + ^ — -jp dr, 


P(n) dr(n) 


(4.51) 


The functions P(n) and Q( tj) are defined In [39] and details on the 
derivation are available in [40]. The Integrals are reduced to finite 
sums by discretizing r(y). Thus, Eq. (4.51) reduces to the following: 


1 N , P( V , 

e i<*i-yf> 0 > -Ti - -y^l (4 - s2 > 


By writing ATfr^) 3 r(y k+1 ) “ Tty^), ^ 4,S2 ^ be expressed as 

N 

e,(x,,y f ,0) 3 Z r(y.)A(i,k) (4.53) 

111 k=l K 

This is the required expression for the downwash due to the r(y) 
distribution. A similar expression is obtained from the second Integral 
on the right-hand side of Eq. (4.43), i.e.. 
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N 

^E^jty^Afi.k) (4.54) 

The vortex sheet behind the wing is assumed to extend downstream from 
the wing trailing edge. The continuous distribution y(x,y) Is replaced 
with a finite number (M) of discrete vortices of strength y^y) 
located at x^, Jt=*l,2,...,M behind the wing. Each of these vortices is 
further represented In the same manner as r(y) Is discretized, so that 
the jet sheet behind the wing Is represented essentially by a vortex 
lattice. The downwash computed from this system of vortices is 

expressed as 

M N 

Mx,,y.,0) -EE Y,(yJB(1,j,lc,A) (4.55) 

J 1 J 1-1 k«l * K 

Equations (4.53) - (4.55) together provide an expression for the 
downwash angle in terms of the unknown r(y k ) , rj(y k ) and y^fy^), 

k=l,2,...,N; 1=1,2, ... ,M. This expression is used in the discretized 
versions of Eqs. (4.44) - (4.46), and the resulting set of simultaneous 
equations are solved for the unknowns. • • • 

The lift-curve slope required in Eq. (4.44) is to be obtained by 
solving Eq. (4.18) or the equivalent Eqs. (4.12), (4.13a) and (4.13c). 

However, in the present analysis, a potential flow method described in 

Chap. 2 has been employed for this purpose. As noted earlier, this 
method has been found to give results which compare reasonably well with 
the solution of Euler equations. 

The analysis presented thus far is applicable to incompressible 
flows. A limited extension is made to account for the compressibility 

7 ~ 

effects by introducing the Prandtl-Glauert factor p = /(1-MJ at 
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appropriate places In the expressions for «, Eq. (4.43), as well as in 
the lift-curve slope in Eq. (4.44). A rigorous treatment of the 
nonuniform compressible flow past a wing section would require the 
solution of the Euler equations coupled with the energy equation. 

However, it is observed that the nonuniformities in the slipstreams at 

high speeds would be relatively smaller than at low speeds. Hence, it 
is hoped that the Prandtl-Glauert correction, which is known to provide 
excellent corrections to the lift of airfoils in uniform subcritical 
flow, will provide a reasonable correction in the present case where the 
approach flow is not uniform. 

4.6 Results And Discussion 

The spanwise lift distribution and total lift have been computed 
for several configurations and compared with available experimental 
data. As the first example, the configuration tested by Stuper [41] is 
considered. In these experiments, the wing was rectangular and spanned 
the wind tunnel walls. The ratio of chord to geometric span was 5.25. 
A slipstream simulator, was placed centrally ahead of the wing. The 

ratio of the slipstream radius to the wing chord was 0.3, and the 
velocity ratio of the slipstream was 1.36. Figure 4.3 shows the 

spanwise distribution of the incremental lift non-dimensionallzed with 
the incremental lift as obtained by the strip theory. Present 

theoretical results are also shown in the figure. It is observed that 
there is a good agreement between the two results, particularly in the 
.region, of the slipstream. Viscous effects at the jet boundary migh^ 
have contributed to the discrepancy around (y/r)=1.0. Experimentally 
measured large values of incremental lift for (y/r)>1.0 could not be 


Ai/(A£) 




(y/r) 


Fig. 4.3 Additional lift distribution due to slipstream on a 
two-dimensional wing at a « 12 degrees. 
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explained, but are suspected to be partly due to the tunnel wall 
Interference. 

As the second example, the configuration tested by Possio [42] is 
considered. In this example, the wing had an aspect ratio of 6.5 and a 
taper ratio of 0.5. A single model propeller was placed ahead of the 
wing. Several sets of data are available with varying velocity ratio 
and angle of attack. Figure 4.4 shows the spanwise lift distribution 
for the wing alone, and for the wing with the slipstreams of velocity 
ratio 1.4 and 1.6. In these tests only the Incremental total lift 
coefficients were measured. The present theoretical values of 
AC l = 0.068 and 0.098 compare very well with the corresponding 
experimental values of 0.07 and 0.10. 

In the next example considered, the wing (tested by Maarsingh [17]) 
had a rectangular planform with an aspect ratio of 7.9. Two simulators 
were placed on each side of the wing at 42 percent and 79 percent of the 
semispan. The simulated slipstream radius was 0.28 times the wing 
chord. -Test results are available for the wing alone and for the wing 
with slipstreams having velocity ratios of 2.0 and 3.0. Figure 4.5 
shows the spanwise lift distribution for the test configurations 
obtained by the present analysis. The spanwise distribution of the 
incremental lift (Ac ^ due to the jets for the two velocity ratios, 
compared with the experimental values in Fig. 4.6. The results computed 
by Maarsingh [16] using Ting's method [15] are also shown in this 
figure. It should be noted that the present method degenerates into 
Ting's, method for unswept wings. However, there are some minor 
differences between the present results and those of Maarsingh [16], 
These differences are suspected to be primarily due to the Inaccurate 
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Fig. 4.4 Spanwlse lift distribution with and without the 
slipstream on a tapered wing of AR a 6.5, 

X = 0.5 at a = 9 degrees. 
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Fig, 4.5 Spanwise lift distribution with and without the 
. slipstream on a rectangular wing of AR = 7.9 at 


« = 9,9 degrees. 
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Fig. 4.6 Spanwise distribution of Incremental lift due to the 
slipstreams on a rectangular wing of AR = 7.9 at 
a = 9.9 degrees. 
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Fig. 4.9 Spanwise load distribution on the swept wing at 
j * 0.8 and a =* 3.28 degrees. 
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outboard stations; but the agreement is not as good at stations close to 
the wing root. It should be noted that the test configuration had a 
body at the center, and the wing planform was modified near the wing 
root. These details were not simulated in the present analysis. Some 
of the differences between the present results and the experimental data 
at the inboard stations may be attributed to these modeling 
difficiencies. 

The effect of the jet without swirl on the spanwise load 
distribution at Mj* 0.8 and a <* 3 degrees is illustrated in Fig. 4.10. 
The experimental data, the results obtained by the present analysis, and 
the solution of Euler equations [26] are Included In this figure. It 

should be noted that the jet is a rather weak jet with a maximum total 

pressure ratio of 1.075. As a result, it does not modify the load 
distribution significantly. It can be observed from the figure that the 
Incremental lift values obtained by the present analysis agree closely 
with the experimental data. A comparison of the Euler solution [27] 

with experimental data shews a trend similar to that which was observed 
earlier - at the wing Inboard stations there is good agreement, but at 
the outboard stations there is some discrepancy. Nevertheless, the 
incremental lift values obtained by this method also compare well with 
the experimental data. 

The effect of the swirl in the jet on the spanwise load 

distribution at M m = 0.8 and a * 3 degrees is illustrated in Fig. 4,11. 
The experimental data are available for the wing with jet having nominal 
swirl angles of 7 degrees up inboard, and 7 degrees down inboard. These 
experimental data are shown In Fig. 4.11a and 4.11b for the up inboard 
and down inboard swirl distributions respectively. Present results and 
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Fig. 4.10 Spanwise load distribution on the swept wing with and 
without the simulator at 0.8 and a» 3 degrees. 
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Fig. 4.11 Effect oi the swirl on the spanwise load distribution 
on the swept wing at M.. = 0.8 and a = 3 degrees. 
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the computed results frc.i; [ 26 ] are also Included in these figures which 
bring out the interesting effect of the swirl in the slipstream on the 
spanwise load distribution. Quantitative agreement between the Euler 
solution [ 26 ] and the experiments is good at stations close to the wing 
root; the agreement is not as good at the outboard stations. In 
contrast, a comparison of the present results with the experiments shows 
a good agreement at the outboard stations but the agreement is not as 
good at stations close to the wing root. This discrepancy in the 
results near the wing root may be attributed to the differences between 
the theoretical model used in the present analysis and the actual wind 
tunnel model. A closer examination of figures reveals that the present 
results agree with experimental data qualitatively; but quantitatively, 
the present analysis seems to slightly overpredict the effect of the 
swirl. 



Chapter 5 
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CONCLUDING REMARKS 


The small perturbation potential flow theory has been applied to 
the problem of determining the chordwlse pressure distribution and lift 
of thin airfoils in nonuniform parallel streams with stepped velocity 
profiles. The method has been extended to the case of an undisturbed 
stream having a given smooth velocity profile with no velocity 
discontinuities. The analysis is based on the method of images, and 
allows for potential disturbances in a rotational undisturbed flow. 
Several examples are considered and the present potential solutions 
obtained by different approximations are compared with the Euler 
solutions. The results Indicate that although approximate, the present 
method yields results which bring out the effect of the interference 
problem, while avoiding the need to solve the Euler equation. In the 
present analysis the thickness effects have been neglected for the 
purpose of computing the airfoil pressure distribution. The effect of 
airfoil thickness can In principle be Included, but would require 
further analysis. 


The classical lifting-line theory for the interference of wing and 
propeller slipstream has been modified in Chap. 3. The classical theory 
assumes that the slipstream is in the form of a circular cylinderical 
jet with a velocity jump across this cylinder. In the present analysis, 
this assumption has been dropped, and the velocity distribution is 


.r-; 



assumed to have a given smooth profile with no velocity 
discontinuities. The lift distribution In the examples considered here 
demonstrated that whereas the total lift may not be affected by the 
assumption of a uniform velocity profile for the slipstream, the 
distributions could be very different in the two cases. The present 
modification to the classical theory brings out the effects of the 
nonuniformity on the spanwise lift and indeed drag distributions. 

The method of analysis developed here can also be extended to the 
case of a wing in multiple nonoverlapping slipstreams. However, It is 
well known that the lifting line theory overpredicts the lift due to the 
slipstream. This is due to the fact that, while accounting for the 
effects of the slipstream boundary, the trailing vortices are assumed to 
extend to Infinity in both directions. This discrepancy was recognized 
and corrected In the development of lifting surface theories. The 

lifting surface theories produce good results, but are applicable only 

to slipstream in which the velocity is uniform. The present method can 
be extended to the lifting surface theories as well as to the lattice 
methods applied to the wing-slipstream interference problem. 

The problem of the interference of multiple propeller slipstreams 
with large aspect ratio swept wings has been treated in Chap. 4 by the 
method of asymptotic expansions. Although this is only a first order 
theory, the important second order contributions from the propeller 
slipstream have been included in the analysis. The method is quite 
versatile in the sense that it can handle slipstreams of any given 

velocity distribution. The time dependence of the velocity field in the 
slipstream is neglected. The increased axial velocity and the swirl 
component in the slipstream are of primary concern in the present 
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analysis. One of the assumptions that the spanwise extent of the 
slipstream was of the order of the wing span, was not satisfied in the 
examples considered; but the results do not seem to be affected 
noticeably. Viscous effects are neglected; but the compressibility 

effect is accounted for by the Prandtl-Glauert factor. 

Several examples are considered for which experimenal data are 

available, and in each case the spanwise lift distribution is obtained 
by the present method, and compared with the corresponding experimental 
data. In most of the cases the agreement between the two results is 
very good. In the first example, the differences between the 
experimental data and the present results are suspected to be partly due 
to the tunnel wall interference in the measurements. In the last 
example considered, the experimental configuration had a body on the 
centerline; also the wing planform was modified near the wing root. 

These details were not simulated in the present analysis. It is 
suspected that these differences have resulted in the difference in the 
two sets of results. 

In the present analysis, the effect of the swirl in the slipstream 
is accounted for' by assuming that the swirl only changes the local angle 
of attack. This approach yields results which compare well 
qualitatively with the experimental data. However, the method seems to 
slightly overpredict the effect of swirl on the spanwise load 

distribution. 

A comparison was made between the present results and the numerical 
solution of Euler equations. The numerical solution of Euler equations 
requires a large computer memory and extensive computing effort. For 



example, the computer program developed by Whitfield and Jameson [24] 
required 900,000 words of memory, and 341 seconds computing time on 
Cray-15 computer for a relatively coarse (96x16x16) grid. This grid 
provided results for only 6 stations along the wing semispan. This 
method, however, provides details of the flowfield in the entire 
computational domain. In contrast, the present method Is simple, and 
requires relatively small computational effort. For example, the 
computer code developed to compute the spanwlse load distribution using 
the present method required 130,000 words of memory, and about 40 
seconds of computing time on CDC Cyber 175 computer, and provided the 
results at 40 stations along the wing semispan. The spanwise lift 

distribution obtained by the present method compared favorably with 
experimental data. 
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